o 

(N 



o 

e 

•4-* 



o 
o 



> 
o 
q 

On 
O 
O 



Impact of degree heterogeneity on the behavior of trapping in Koch networks 
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Previous work shows that the mean first-passage time (MFPT) for random walks to a given hub node (node 
with maximum degree) in uncorrelated random scale-free networks is closely related to the exponent 7 of 
power-law degree distribution P(k) ~ fc -7 , which describes the extent of heterogeneity of scale-free network 
structure. However, extensive empirical research indicates that real networked systems also display ubiquitous 
degree correlations. In this paper, we address the trapping issue on the Koch networks, which is a special 
random walk with one trap fixed at a hub node. The Koch networks are power-law with the characteristic 
exponent 7 in the range between 2 and 3, they are either assortative or disassortative. We calculate exactly 
the MFPT that is the average of first-passage time from all other nodes to the trap. The obtained explicit 
solution shows that in large networks the MFPT varies lineally with node number A, which is obviously 
independent of 7 and is sharp contrast to the scaling behavior of MFPT observed for uncorrelated random 
scale-free networks, where 7 influences qualitatively the MFPT of trapping problem. 
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As a fundamental dynamical process, random 
walks have received considerable interest from 
the scientific community. Recent work shows 
that the key quantity — mean first-passage time 
(MFPT) for random walks to a given hub node 
(node with highest degree) on uncorrelated ran- 
dom scale-free networks is qualitatively reliant on 
the heterogeneity of network structure. However, 
in addition to the power-law behavior, most real 
systems are also characterized by degree correla- 
tions. In this paper, we study random walks on a 
family of recently proposed networks — Koch net- 
works that are transformed from the well-known 
Koch curves and have simultaneously power-law 
degree distribution and degree correlations with 
the power exponent of degree distribution lying 
between 2 and 3. We explicitly determine the 
MFPT, i.e., the average of first-passage time to a 
target hub node averaged over all possible start- 
ing positions, and show that the MFPT varies 
linearly with node number, independent of the in- 
homogeneity of network structure. Our result in- 
dicates that the heterogeneous structure of Koch 
networks has little impact on the scaling of MFPT 
in the network family, which is in contrast with 
result of MFPT previously reported for uncorre- 
lated stochastic scale-free graphs. 



a ' Electronic mail: zhangzz@fudan.edu.cn 



I. INTRODUCTION 



In the past decade, a lot of endeavors have been de- 
voted to characterize the structure of real systems from 
the view point of complex networks^— , where nodes rep- 
resent system elements and edges interactions or relations 
between them. One of the most important findings of 
extensive empirical studies is that a wide variety of real 
networked systems exhibit scale-free behavior—, charac- 
terized by a power-law degree distribution P(k) ~ k 1 
with degree exponent 7 lying in the interval of [2,3]. Net- 
works with such broad tail distribution are called scale- 
free networks, which display inhomogeneous structure en- 
coded in the exponent 7: the less the exponent 7, the 
stronger the inhomogeneity of the network structure, and 
vice versa. The heterogeneous structure critically influ- 
ences many other topological properties. For instance, it 
has been shown that in uncorrelated random scale-free 
networks with node number N (often called network or- 
der), their average path length (APL)£ relies on 7^: for 
7 = 3, d(N) ~ In A; while for 2 < 7 < 3, d(N) ~ In In A. 

The power-law degree distribution also radically af- 
fects the dynamical processes running on scale-free net- 
works^, such as disease spreading^, percolation^, and 
so on. Amongst various dynamics, random walks are an 
important one that have a wide range of applications^ - — 
and have received considerable attention^—. Recently, 
mean first-passage time (MFPT) for random walks to 
a given target point in graphs, averaged over all source 
points has been extensively studied^— . A striking find- 
ing is that MFPT to a hub node (node with highest de- 
gree) in scale-free networks scales sublinearly with net- 
work order— ~—, the root of which is assumed to be the 
structure heterogeneity of the networks. In particularly, 
it has been reported 28 that in uncorrelated random scale- 




FIG. 1. (Color online) Iterative construction method for the 
Koch networks. 



free networks, the MFPT F(N) scales with network or- 

7-2 

der N as F(N) ~ jVt- 1 . However, real networks exhibit 
ubiquitous degree correlations among nodes, they are ei- 
ther assortative or disassortative^. Then, an interesting 
question arises naturally: whether the relation governing 
MFPT and degree exponent 7 in uncorrelated scale-free 
networks is also valid for their correlated counterparts. 

In this paper, we study analytically random walks in 
the Koch network o 30 ' 31 that are controlled by a positive- 
integer parameter to. This family of networks is scale- free 
with the degree exponent 7 lying between 2 and 3, and 
it may be either disassortative (m > 1) or uncorrelated 
(to = 1). We focus on the trapping problem, a particular 
case of random walks with a fixed trap located at a hub 
node. We derive exactly the MFPT that is the average 
of first-passage time (FPT) from all starting nodes to 
the trap. The obtained explicit formula displays that in 
large networks with N nodes, the MFPT grows linearly 
with N, which is independent of 7 and, showing that 
the structure inhomogeneity has no quantitative influ- 
ence on the MFPT to the hub in Koch networks, which 
lies in their symmetric structure and other special fea- 
tures and is quite different from the result previously re- 
ported for uncorrelated random scale-free networks. Our 
work deepens the understanding of random walks occur- 
ring on scale-free networks. 



II. CONSTRUCTION AND PROPERTIES OF KOCH 
NETWORKS 

The Koch networks governed by a parameter to are de- 
rived from the famous Koch curve o 32 i 33 , which are con- 
structed in an iterative wa y 30 ' 31 . Let K m _ t denote the 
Koch networks after t iterations. Then, the networks can 
be generated as follows: Initially (t = 0), K m fi is a trian- 
gle. For t > 1, K mA is obtained from K m> t-i by adding 
to groups of nodes for each of the three nodes of every 
existing triangles in K m ^-i- Each node group consists 
of two nodes, both of which and their "father" node are 
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FIG. 2. (Color online) A network corresponding to the case 
of m = 2. 

connected to one another shaping a new triangle. That 
is to say, to get K m t from K m< t-i, one can replace each 
existing triangle in K m ^-i by the connected clusters on 
the right-hand side of Fig. [TJ Figure [5] show a network 
corresponding to to = 2 after several iterations. 

By construction, the total number of triangles L/\{t) 
at iteration t is L&(t) = (3m + 1)', and the number of 
nodes created at iteration t is L v (t) — 6mL&(t — 1) = 
6to (3to + 1)* — 1 . Then, the total number of nodes N t 
present at step t is 

t 

JV t =5^£ w (ti) = 2(3m + l)* + l. (1) 
t«=o 

Let ki(t) be the degree of a node i at time t, which 
is added to the networks at iteration (step) ij (U > 0). 
Then, fcj(tj) = 2. Let L/\(i,t) denote the number of 
triangles involving node i at step t. According to the 
construction algorithm, each triangle involving node i at 
a given step will give birth to to new triangles passing by 
node i at next step. Thus, L/\(i, t) = (to + 1) L&(i, t — 
1) = (to + l)* - **. Moreover, it is easy to have fcj(t) = 
2L A (i,t), i.e., 

k l {t)=2L A {i,t)=2{m+l) t - t \ (2) 

which implies 

h[t) = {m + l)k i {t-l), (3) 

The Koch networks present some common features of 
real systems-^. They are scale-free, having a power-law 
degree distribution P{k) ~ k^ 1 with 7 = 1 + ^(m+i) 
longing to the range between 2 and 3. Thus, parameter 
to controls the extent of heterogeneous structure of Koch 
networks with larger to corresponding to more heteroge- 
neous structure. They have small- world effect with a low 
APL and a high clustering coefficient. In addition, their 
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FIG. 3. (Color online) labels of all nodes in 7^2,2- 



degree correlations can be also determined. For m = 1, 
they are completely uncorrelated; while for other values 
of m, the Koch networks are disassortative. 



III. RANDOM WALKS WITH A TRAP FIXED ON A 
HUB NODE 

After introducing the construction and structural 
properties of the Koch networks, we continue to inves- 
tigate random walks 34 performing on them. Our aim 
is to uncover how topological features, especially degree 
correlations influence the behavior of a simple random 
walk on Koch networks with a single trap or a perfect 
absorber stationed at a given node with highest degree. 
At each step the walker, located on a given node, moves 
uniformly to any of its nearest neighbors. To facilitate 
the description, we label all the nodes in K m t as fol- 
lows. The initial three nodes in K m have label 1, 2, 
and 3, respectively. In each new generation, only the 
newly created nodes are labeled, while all old nodes hold 
the labels unchanged. That is to say, the new nodes are 
labeled consecutively as M+l, M + 2, ...,M + AM, with 
M being the number of all pre-existing nodes and AM 
the number of newly created nodes. Eventually, every 
node has a unique labeling: at time t all nodes are la- 
beled continuously from 1 to N t = 2(3™ + 1)* + 1, see 
Fig. [3] We locate the trap at node 1, denoted by It- 

We will show that the particular selection of the trap 
location makes it possible to compute analytically the 
relevant quantity of the trapping process, i.e., mean first- 
passage time. Let if denote the first-passage time of 
node i in K m ^ except the trap ir, which is the expected 
time for a walker starting from i to visit the trap for the 

first time. The mean of FPT if over all non-trap nodes 
in K m t is MFPT, denoted by (F) t , the determination of 
which is a main object of the section. To this end, we 
first establish the scaling relation governing the evolution 
of if with generation t. 



TABLE I. First-passage time if for a random walker start- 
ing from node i in K m ,t for different t. Note that thanks to 
the symmetry, nodes in the same column are equivalent to 
one another, since they have the same FPT. 



t/i 


2-3 


4-7 8-15 


16 - 27 28 - 51 


52-67 


68 - 99 





2 














1 


14 


2 


16 










2 


98 


14 


112 


2 


100 


16 


114 


3 


686 


98 


784 


14 


700 


112 


798 


4 


4802 


686 


5488 


98 


4900 


784 


5586 


5 


33614 


4602 


38416 


686 


34300 


5488 


39102 



A. Evolution scaling for first-passage time 

We begin by recording the numerical values of if for 
the case of m = 2. Clearly, for all t > 0, if) = 0; for 
t = 0, it is trivial, and we have if } = if) = 2. For 

t > 1 , the values of if can be obtained numerically but 
exactly via computing the inversion of a matrix, which 
will be discussed in the following text. Here we only give 
the values of computation. In the generation n = 1, by 
symmetry we have F 2 (1) = F 3 (1) = 14, if) = if ] = 
if) = if) = 2, and if) = if) = . . . = F$ = 
16. Analogously, for t = 2, the numerical solutions are 
if) = if) = 98, if' = if) = if) = if) = 14, 
F 8 (2) = if' = . . . = if) = 112, if) = if 2) =...= 

77.(2) _ 9 p (2) _ p (2) _ _ p (2) _ 1m p (2) _ p (2) 
r 27 ~ z ' r 2S ~ r 29 — • • • — ^51 — J-UU, r 52 — ^53 — 

. . . = F 6 (2) = 16, and F 6 (2) = F 6 (2) = = F 9 (2) = 114. 

Table Q] lists the numerical values of if for some nodes 
up to t = 5. 

The numerical values reported in Table U show that 
for any node i, its FPT satisfies the relation if +1 ) = 

(3m + l)if \ In other words, upon growth of Koch 
networks from generation t to t + 1, the FPT of any 

(5) 

node increases to 3m + 1 times. For example, F 2 = 
7if ) = 7 2 if ' = 7 3 if ) = 7 4 if ) = 7 5 if ) = 33614, 
F 8 (5) - 7if ) = 7 2 F 8 (3) = 7 3 F 8 (2) = 7 4 if) = 38416, and 
so forth. This scaling is a basic property of random walks 
on the family of Koch networks, which can be established 
based on the following arguments. 

Examine an arbitrary node i in the Koch networks 
K m j- Equation shows that upon growth of the net- 
works from generation t to t + 1, the degree ki of node i 
grows by m times, i.e., it increases from ki to (m + l)ki. 
Let A denote the FPT for going from node i to any of 
its ki old neighbors, and let B be FPT for starting from 
any of the mki new neighbors of node i to one of its fc; 
old neighboring nodes. Then the following equations can 
be established (see Fig. 2]): 
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FIG. 4. (Color online) Growth of first-passage time in going 
from K m ,t to Km : t+i in the case of m — 2. Node i £ K m , t has 
ki neighbor nodes in generation t (O) an d mki new neighbor 
nodes in generation t + 1 (□). A new neighbor of node i has 
a degree of 2, and is simultaneously linked to another new 
neighbor of i. 




FIG. 5. (Color online) Illustration showing the relation of the 
first-passage times for each pair of two new nodes (v x and w x 
with x=l,2,. . . , or m) and the old node u as one point of the 
triangle generating the new nodes. 



which yield A = 3m + 1. This indicates when the net- 
works grow from generation t to t + 1, the FPT from any 
node Fm,t) to any node j (j £ K 7riit +i) increases on 



average 3m times. Then, we have F t 



(t+i) 



(3m + l)F t 



(*) 



For explanation, see Refs . 35 ' 36 and related references 
therein. The obtained relation for FPT is very useful 
for the following derivation of MFPT. 



B. Explicit expression for mean first-passage time 

Having obtained the scaling dominating the evolution 
for FPT, we now draw upon this relation to determine 
the MFPT, with an aim to derive an explicit solution. 
For the sake of convenient description of computation, 
we represent the set of nodes in K mi t as 0t, and denote 
the set of nodes created at generation t by 6 t . Evidently, 
the relation Q t = Q t U Qt-i holds. In addition, for any 
r < t, we dchnc the two following variables: 



F. 



(t) 



E 



F 



(*) 



and 



Then, we have 



F, 



E ^ (t) - 



(5) 



(6) 



F 



(t) 



F 



it) 



(Zm + VFtxlt + FXot, (7) 



t,tot t— l,tot i,tot 

and 



(*) 



(F)t 



F 



(t) 



AT* — 1 



Thus, to explicitly determine the quantity {F) t , one 
should first find F^ ot , which can be reduced to deter- 
mining F t fj ot . Next, will show how to solve the quantity 



F 



U) 



By construction, at a given generation, for each trian- 
gle passing by node w, it will generate m new triangles 
involving u (see Fig. [5]). For each of the m new triangles, 
the first-passage times for its two new nodes (v x and w x ) 
and that of its old node u follow the relations: 



F{v x ) = 1 + 4 [F(0 + F(u)] , 
F(w x ) = l + \[F(v x ) + F(u)} . 



(9) 



In Eq. ©, F(s) represents the expected time of a particle 
to first visit the trap node, given that it starts from node 
s. Equation ([9]) yields 



F(v x ) + F{w x ) = A + 2F(u). 



(10) 



Summing Eq. ([TO]) over all the L A (t) = (3m + 1)' old 
triangles pre-existing at the generation t and the three 
old nodes of each of the L/s(t) triangles, we obtain 

F t ( *+Jot = 3 • 4 • mL A (t) + J2 (2mL A (i,i) • if 

iee t 

= 12m(3m + 1)* + 2mF t % 1) + 2m(m + l^l+Jot 
+ . . . + 2m(m + l) 4 " 1 ^ + 2m(m + l)*^^ . 

(11) 

For instance, in K 2 ,2 (see Fig. [3]), tot can be expressed 
as 



99 

41 = E^ (2) = 1176 + 12 A (2) 

2 = 16 



„., . 36F$ ot . (12) 



Now, we can determine F^ ot through a recurrence re- 
fg\ lation, which can be obtained easily. From Eq. (fl2"|) . it is 



not difficult to write out F^Jtot as 



F 



(t+2) 



12m(3m + l) t+i + 2mF t ( ^ 2 t ) ot + 2m(m + l)F t 
+ ... + 2m(m + l)*^^ + 2m(m + l) t+1 F ( * 



t+2, tot 



(13) 
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Equation (IT3l) minus Eq. (jlll) times (m + l)(3m+l) and 
using the relation F 7 



0+2) 



(3m+ we have 



F 



(t+2) 



t+2,tot 



(3m + 1)^,1 



12m 2 (3m + 1 
(i) 



,t+i 



(14) 



Making use of the initial condition F{ tot = 24m + 20m, 
Eq. (fT4|) is solved inductively to yield 

F t ( * } ot = 4m(3m + l)'" 1 + (24m 2 + 16m)(3m + 1) 2 *~ 2 . 

(15) 

Inserting Eq. (fTo*)) for F^ ot into Eq. flSJ), we have 



i% t = (3m + l)F^_tL 
4m(3m + 1)'" 



(24m 2 + 16m)(3m + l) 2t ^ 2 . 

(16) 



Since -Fotot = 4, we can resolve Eq. (fl~6|) by induction to 
obtain 

^t.tot = ^(3m+l) t " 1 [(6m + 4)(3m + l)* + 3mi + 3m-l] . 

(17) 

By plugging Eq. (|17|) into Eq. we obtain the closed- 
form solution to the MFPT for random walks on the 
Koch networks with an immobile trap stationed at a hub 
node: 



(F)t 



3(3m+ 1) 



[(6m + 4) (3m + 1)* + 3mi + 3m - 1] . 

(18) 
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FIG. 6. (Color online) Mean first-passage time {F)t as a func- 
tion of the generation t on a semilogarithmic scale for different 
values of m. The empty symbols represent the numerical re- 
sults obtained by direct calculation from Eq. (|20p ; while the 
filled symbols correspond to the rigorous values provided by 
Eq. US]). 



where (t) is the expected times that the walk visit 
node j, given that it starts from node *22. Using Eq. (fH)]) 

we can determine F^ numerically but exactly for differ- 
ent non-trap nodes at various generation t, as listed in 
Table D 

By definition, the MFPT (F) t that is the mean of F- (t) 
over all initial non-trap nodes in K m j reads as 



C. Numerical calculations 

We have corroborated our analytical formula for 
MFPT provided by Eq. (fT5)l against direct numerical cal- 
culations via inverting a matri:»2I Indeed, the Koch net- 
work family K m> t can be represented by its adjacency 
matrix A t of an order N t x N t , the element ay(t) of 
which is either 1 or defined as follows: Oij(i) = 1 
if nodes i and j are directly connected by a link, and 
a ij(t) — otherwise. Then the degree, di(t), of node 
i in K mt t is given by d t (t) = £V ' ^ji 1 )-, tne diag- 
onal degree matrix Z t associated with K m:t is Z 4 = 
diag(di(t), d<i{t), . . . , di(t), . . . , dpj t (t)), and the normal- 
ized Laplacian matrix of K m j is provided by L t = 
l t — Z t " 1 A t , in which I t is the Nt x N t identity matrix. 

Note that the random walks considered above is in 
fact a Markovian process, and the fundamental matrix of 
Markov chain representing such unbiased random walks 
is the inverse of a submatrix of L t , denoted by L t that 
is obtained by removing the first row and column of L t 
corresponding to the trap node. According to previous 
result^, the FPT F^ can be expressed by in terms of 
the entry l^- (i) of L^ 1 as 



F 



N t 
J=2 



(19) 



1 Nt 

1 i=2 



(t) 



Nt N 



2(3m 



(20) 



i=2 j=2 



In Fig. [5J we compare the analytical results given by 
Eq. (fT8"]) and the numerical results obtained by Eq. (|2H)) 
for various i and m. Figure [6] shows that the analyti- 
cal and numerical values for (F) t are in full agreement 
with each other. This agreement serves as a test of our 
analytical formula. 



D. Dependence of mean first-passage time on network 
order 

Below we will show how to express {F) t as a function of 
network order Nt, with the aim of obtaining the relation 
between these two quantities. Recalling Eq. ([T]), we have 

(3m + and t = ln( g 3 -^- 1 ) n2 - Thus > E ^ dUD 

can be recast as in terms of Nt as 



_ 2(3m + 2) 
{F)t ~ 3(3m + l) (7Vt ~ 1} 
2(3m- 1) 
3(3m + l) ' 



2m[\n(N t - 1) - In 2] 
(3m + 1) ln(3m + 1) 



(21) 
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In the thermodynamic limit (N t — > oo), we have 



(22) 



showing that the MFPT grows linearly with increasing 
order of the Koch networks. Equations (|21|) and (|22|) 
imply that although for different m the MFPT of whole 
family of Koch networks is quantitatively different, it ex- 
hibits the same scaling behavior despite the distinct ex- 
tent of structure inhomogeneity of the networks, which 
may be attributed to the symmetry and particular prop- 
erties of the networks studied. 

It is known that the exponent 7 characterizing the in- 
homogeneity of networks affects qualitatively the scaling 
of MFPT for diffusion in random uncorrelated scale-free 
networks^,. Concretely, in random uncorrelated scale- 
free networks with large order N, the MFPT F(N) grows 

sublinearly or linearly with network order as F(N) ~ 

-7-2 

N-1- 1 for all 7 > 2, which strongly depends on 7. How- 
ever, as shown above, in the whole family of Koch net- 
works, the MFPT displays a linear dependence on net- 
work order, which is independent of 7, showing that the 
inhomogeneity of structure has no quantitative impact 
on the scaling behavior of MFPT for trapping process 
in Koch networks. Our obtained result means that the 
scaling observed for MFPT in the literature 28 is not a 
generic feature of all scale-free networks, at least it is not 
valid for the Koch networks, even for the case of m = 1 
when network is uncorrelated. 



IV. CONCLUSIONS 

Power-law degree distribution and degree correlations 
play a significant role in the collective dynamical behav- 
iors on scale-free networks. In this paper, we have inves- 
tigated the trapping issue, concentrating on a particular 
case with the trap fixed on a node with highest degree on 
the Koch networks that display synchronously a heavy- 
tailed degree distribution with general exponent 76 [2,3] 
and degree correlations. We obtained explicitly the for- 
mula for MFPT to the trapping node, which scales lin- 
eally with network order, independent of the exponent 
7. Our result shows that structural inhomogeneity of the 
Koch networks has no essential effect on the scaling of 
MFPT for the trapping issue, which departs a little from 
that one expects and is as compared with the scaling 
behavior reported for stochastic uncorrelated scale-free 
networks. Thus, caution must be taken when making a 
general statement about the dependence of MFPT for 
trapping issue on the inhomogeneous structure of scale- 
free networks. Finally, it should be also mentioned that 
both random uncorrelated networks and the Koch net- 
works addressed here cannot well describe real systems, 
future work should focus on trapping problem on those 
networks better mimicking realties. Anyway, our work 
provides some insight to better understand the trapping 
process in scale-free graphs. 
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